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Abstract 
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We consider the conversions of active to sterile Majorana neutrinos v a and 

ID ■ 

. u s , due to neutrino transition magnetic moments in the presence of random 

G\ 

o 

a simple Schrodinger-type evolution equation, we derive a stringent constraint 

ON ' 

on the corresponding transition magnetic moments and display it as a function 
Qh| of the domain size and field geometry. For typical parameter choices one gets 

Q_i' limits much stronger than usually derived from stellar energy loss considera- 

e-j . tions. These bounds are consistent with the hypothesis of seeding of galactic 



magnetic fields (r.m.f.) generated at the electroweak phase transition. From 



magnetic fields by primordial fields surviving past the re-combination epoch. 
^ . We also obtain a bound on active-sterile neutrino transition magnetic moments 

from supernova energy loss arguments. For r.m.f. strengths in the range 10 7 
to 10 12 Gauss we obtain limits varying from LL v as < 10 13 /xs to fj,^ s < 10 18 /x#, 
again much stronger than in the case without magnetic fields. 
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1. Introduction 



Recently there has been a renewed interest concerning neutrino propagation in me- 
dia with random magnetic fields, both from the point of view of the early universe 
hot plasma as well as astrophysics 0. It has been shown that random magnetic 
fields can strongly influence neutrino conversion rates and this could have important 
cosmological and astrophysical implications, especially in the case of conversions in- 
volving a light sterile neutrino v s . Indeed, note that the present hints from solar and 
atmospheric neutrino observations 0|| as well as from the COBE data on cosmic 
background temperature anisotropics on large scales HH indicate that a light sterile 
neutrino v s might exist in nature. 

So far the most stringent constraints for the neutrino mass matrix including a 
sterile neutrino species, u s , are obtained from the nucleosynthesis bound on the max- 
imum number of extra neutrino species that can reach thermal equilibrium before 
nucleosynthesis and change the primordially produced helium abundance ||. This 
has been widely discussed in the case of vanishing magnetic field || and for the 
case of a large random magnetic field that could arise from the electroweak phase 



transition and act as seed for the galactic magnetic fields [ |10] . 

The effect of active-sterile neutrino conversions in a supernova has also been dis- 



cussed, both in the case where no magnetic field is present [Tl[, as well as in the 
presence of random magnetic field as large as 10 16 Gauss fL2|, following a suggestion 



made in ref. [ 13 



In this paper we focus on the important effect of relatively small random magnetic 
fields on the active sterile neutrino conversion rates when nonzero Majorana neutrino 
transition magnetic moments are taken into account |I^JT5[1 . We apply this to the 
case of v a to v s conversions in the early universe hot plasma, as well as in a supernova, 
showing how their effect can place limits that are substantially more stringent than 
those that apply in the absence of a magnetic field. First we focus on the nucleosyn- 
thesis constraint on active to sterile transition magnetic moments, and display this 
constraint as a function of the domain size and field geometry. For typical choices 
we obtain [i v as < 10~ 15 fis, which is much stronger than usually derived from stel- 
lar energy loss considerations neglecting r.m.f. effects. These bounds are consistent 
with the hypothesis of seeding of galactic magnetic fields by primordial fields gener- 
ated at the electroweak phase transition. We also obtain a bound on active-sterile 
neutrino transition magnetic moments from supernova energy loss arguments. For 
modest r.m.f. field strengths in the range 10 7 to 10 12 Gauss we obtain limits that 
vary from jji u as < 10~ 13 yU B to ji u as < lCT 18 /iB, again much stronger than in the case 



without magnetic fields [16 



2. Neutrino conversions 



Let us consider the Schrodinger equation describing a system of two neutrinos, one 
active and one sterile, v a and v s in a plasma with a random magnetic field B(t). In 



general one has a four- dimensional system of equations W% - We will use the ultra- 
relativistic limit and will neglect the corresponding mixing angle, s — 0, c = 1, in 
which case one may write [17] 

±(v a \ = (V a -A + ii eff B [{ (t) fxB ± (t)\/u a \ 
% dt\u s ) \ fjiB ± {t) )\u a )' 

Here V a (for u a , a — e,fi, r) is the active neutrino vector interaction potential. For 
instance, for electron left-handed neutrinos propagating in the early universe hot 
plasma one has 

V = 3.45 x lO-™(-L-y M eV, (2.2) 
while for the case of a supernova one has 

V-4x 10- 6 pi4(3F e + 4Y Ve + 2Y; M - l)MeV. (2.3) 

where pu is the core density in units of 10 14 g/cm 3 . 

In the case of the muon left-handed neutrino the last abundance factor in eq. 
( |2.3|) is changed to f(Y) = Y e — 1 + 4F„ M + 2Y Ue , where we took into account v-v 
forward scattering amplitude in a supernova matter during the main neutrino burst 
time, about 10-20 seconds. 

The term fi e ffB\\ is produced by the mean axial vector current of charged leptons 
in an external magnetic field B\\ = B.q/g. In the case of the early universe hot 
plasma the coefficient /i e // is given as 

T 

while for the case of an ultra-relativistic degenerate electron gas of a supernova one has 

fieff = -8.Qc A x 10- 13 fi B {p Fe /MeV) . (2.5) 

Here ca = =F0.5 is the axial constant (upper sign for electron neutrino, lower one 
for v^ T ). 

Note that even though the quantity /i e // has the dimensions of a magnetic moment 
it is not a real magnetic moment since it is helicity conserving. Indeed, the additional 
energy splitting term obtained in eq. ( |2.1| ) does not lead to any helicity change. In 
contrast, the off-diagonal entries of the Hamiltonian eq. ( |2.1| ) involve the presence 
of the real neutrino transition magnetic moment /i = p u as as well as the transversal 
magnetic field component, B^. 



/i e// = -12c A xlO- 13 ^(— -), (2.4) 



Note also that our initial equation eq. (|2TT| ) is quite general, since it applies to the 
case of Majorana neutrinos. It remains valid also for the description of the active to 
sterile neutrino conversions in the limit where these form a Dirac neutrino, in which 
case A = (m?, — mf) /2q = and [i becomes the usual magnetic moment. 

Denoting fiB^(t) = H±(t) and using the auxiliary functions R = Re((u*u s )) 
and / = Im((u*u s )), we obtain from eq. Q2.1|) the standard system of first order 
differential equations, 

P(t) = -2H ± (t)I(t), 

i(t) = [V - A + Mf B\\{t)]R{t) + H ± (t)(2P(t) - 1), 

R(t) = -[V-A + Mf B {l (t)]I(t), 

from which we derive the integro-differential equation for the neutrino conversion 
probability P„ e ^ Vs (t) = P(t) 

P(t) = -2 jf* cos(^* V{t 2 )dt 2 ) [2P(tx) - l]dtx, (2.6) 

where jj V(t 2 )dt 2 = (V - A)(t - h) + M} Jj B\\{t 2 )dt 2 . 

For r.m.f. domain sizes much smaller than the typical neutrino conversion length 
honv ~ r _1 , i- e - Lq Iconv (see below eq. Q2.11| )) one can average this equation over 
the random magnetic field distribution. The averaged probability (P(t)) = V(U) 
must depend only on even powers of the random field, since ((Bj(t)) 2n+1 ) = 0. 

In what follows we neglect neutrino collisions with dense matter. One can show 
that averaging of the neutrino conversion probability over fast collisions does not 
lead to essential change of our results f\ 

For uncorrelated random magnetic field domains, 

(B l {x)B 3 {y)) = ^5 l3 ^\x-y) ) (2.7) 
we obtain the mean values 

(S||(t)> = (B i .(t)) = <B|,(t)S J .(t 1 ))=0, 

(B ll (t)B ll (t 1 )) = (Bf l )L S(t-t 1 ), (2.8) 

(Bx(t)B±(ti)} = (BDLaSit-h), 

The correlation length A is determined by the domain size L and by the value of 
root mean squared field B = B rms = yj (B 2 ) at the horizon scale L = Ih |T^fl : 

l = n(3-2p) B * ms{lH)L3 °> 

3 Moreover, in [^9| it was shown that in an appropriate kinetic approach for the neutrino spin 
evolution in uncorrelated random magnetic fields the final result does not depend on the collision 
frequency at all 



for p 7^ 3/2, and 

A it L 



for p = 3/2. 

In eq. the mean squared fields {B? { ) = (5 2 )/3, (B\) = 2(B 2 )/3 are given 



by the r.m.s. field B = \J(B 2 }. 

Averaging the master equation eq. ( |2.6| ) over a random magnetic field distribution 
with the use of eq. ( |2.8[ ) we easily obtain the simple differential equation 

V + T ± V = (2.9) 

which has the solution 

V = ^(oo - exp(— ± U)) « (2.10) 
where the magnetic field damping parameter r_i_ is defined as 

r ± = Afi 2 (Bl)Lo = ^ 2 B 2 L . (2.11) 
Note that the result eq. fl2.10| ) obtained from our Schrodinger equation eq. ( p.l|) 



coincides with the analogous ones obtained in [19j for the Dirac neutrino magnetic 
moment, using neutrino spin kinetic equation PD| , or the Redfield equation for the 
density matrix describing a nuclear spin interacting with the lattice vibrations in 
solids [|2T . 



Note also that in the simplified case of negligible neutrino mixing considered here, 
the assumed 5-correlated form of the random magnetic field implies that there is no 
dependence of the averaged neutrino conversion probability on the neutrino energy 
splitting due to the matter term V, to the term p, e ffB\\ or to the mixing parameter A. 



3. Nucleosynthesis constraint 

The rate for producing sterile neutrinos in the early universe hot plasma is given as the 
product of a typical weak interaction rate IV times our averaged conversion probabil- 
ity V . The constraint that follows from nucleosynthesis may be simply estimated as 

T s = T w V<n, (3.1) 

where rv = 4.0G^T 5 is the rate for producing the standard model active neutrinos, 
H = 4.46 x 10~ 22 (T/MeV) 2 MeV is the Hubble parameter, and V is given by eq. 
»• 

Different mechanisms of the magnetic field generation in the early universe give 
the following main phenomenological formula 




where we put the maximal scale is chosen as the horizon scale f] L = ln{t)- 

Note that if we put in eq. ( j3T2"D B ~ T 2 (as in Vachaspati's mechanism, where 
Bo ~ T EW [HD the r - m - s - field i n eq. (|2.10|) depends only on the temperature of the 



universe T (T ~ Tqqd m eq. (|3.3| )) and on two parameters describing the geometry of 
the r.m.f. geometry: (a) the scaling parameter p, and (b) the size of random domain 
L™ n = a(MeV/T). 

Using this equation we get the following constraint | 

< 1.7 x io- 2 y B 

where L n%n is the minimal scale of the random magnetic field domains. Due to the 
copious production of left-handed neutrinos e.g. from pion decays after the quark 
hadron phase transition temperature T ~ Tqqd — 200 MeV this QCD temperature 
plays a crucial role in the estimate of the nucleosynthesis constraint. As a result all 
parameters in eq. (|3.3| ) are evaluated at T ~ Tq C d — 200 MeV. 

The allowed region of transition moments fi u we have obtained from nucleosyn- 
thesis is shown in Fig. 1 as a function of the scale parameter p and domain size Lq. 

There is an independent bound from primordial nucleosynthesis on the strength 
of the random magnetic field. Indeed, it has been shown |23| that a too large r.m.f. 



strength would enhance the rates for the relevant weak processes at Tjvs ~ 0.1 MeV 
resulting in additional helium production. This leads to B < 3 x 10 10 Gauss |]23 
which implies, at the present time flT9| , 

B N s < 1-8 x 2.4 P x 10- 7 - llp Gauss (3.4) 

Clearly, all values of the parameter p between p = (corresponding to uniform field) 
and p = 3/2 (corresponding to 3-dimensional elementary cells) can obey both nucle- 
osynthesis limits on magnetic transition moments eq. ( |3.3|) and on r.m.f. strength, 
eq. (H). 

It is important to stress that, for a wide choice of p and L , our nucleosynthesis 
constraint on the active sterile transition moments is substantially stronger than the 



astrophysical one from supernova 1987A [16 



Until now our arguments are quite general. Now we show that our constraints 
are consistent with the hypothesis of a primordial origin of the observed galactic 
magnetic fields p2|j24| which would require additional restrictions on the parameters 



p and Lo, resulting in a correspondingly stronger limit on the neutrino transition 
magnetic moments. 



i As explained in ref. this follows from the definition of the parameter A for the initial 



5-correlator in eq. (2.7). 

5 Clearly in this approximation eq. (3.1) leads to the same constraint on the transition magnetic 



moment /i — fi a s as obtained previously using kinetic theory for the case of Dirac neutrino diagonal 
magnetic moment fig| . 



This would require L™ m > 10 3 (MeV/T) cm and p < 1. The first restriction 
comes from requiring that primordial fields survive against ohmical dissipation after 
recombination time and could seed the observed galactic magnetic field ||23|| . On 
the other hand, the limit p < 1 follows from eq. ( |3.4|) and the lower dynamo theory 
astrophysical bound on the strength of the galactic magnetic field B seed (T now ) > 1CT 18 
Gauss ||25|| . The corresponding region is hatched in Fig. 1. This corresponds to an 
upper limit on the neutrino transition moment /i < 10~ 16 /zb (left corner of hatched 
region for p = 1). Moreover, in order to be a seed field for the Milky Way or 
Andromeda, one would require B seed (T now ) > 1CT 13 Gauss E5j which corresponds to 



p < 0.5 [EH], leading to a much tighter bound /i < 10 19 /is. 



4. Supernova bounds 

Now let us move to the case of a supernova with magnetic field. If this magnetic 
field is generated after collapse it could be viewed as a random superposition of many 
small dipoles of size L ~ 1 km |TJ|. One can show that this hypothesis is consistent, 
even for large magnetic field strengths, with the non observation of gamma radiation 



from SN 1987A, if the magnetic field diffusion time [13 



t diff ~ l0 2 (B /10 9 G) 2 sec (4.1) 

does not exceed the time necessary for the diffusion of the X-rays through the super- 
nova mantle. This leads to an upper limit on the seed field B , B < 10 12 G. Thus, as 
long as this is fulfilled, the random magnetic field could well influence the SN 1987A 
neutrino burst without being observable at present via X-ray emission []. 

As we will show below, one can derive constraints on the active to sterile Majorana 
neutrino magnetic moments which are so stringent that they become relevant even 
for values of the r.m.s field strength as low as B < 10 9 Gauss which could be quite 
reasonable from the astrophysical point of view. 

First, let us consider the strong small-scale random magnetic fields generated in 
a supernova as in |13|]. For this case we can use the requirement that in the non- 
trapping regime the sterile neutrino can be emitted from anywhere inside the stellar 



core with a rate 12 



^ - 2.8 x 10 55 V PooA S^ < oo, A £ (4.2) 

where pu is the core density in units of 10 14 g/cm 3 and E'ioo is the neutrino energy 

in units of 100 MeV. In the following we will set these parameters to unity. The last 

6 Alternatively, for the case of neutron stars the assumption of random magnetic field domains is 
consistent with their observed magnetic fields in the limit where the small domains merge together 
to form larger ones due to kinematic evolution of the random magnetic field energy transferred down 
the spectra E(k) ~ fc~ 5 / 3 Jl| with the characteristic time t(Lq) ~ L^ 3 



bound corresponds to the maximum observed integrated neutrino luminosity. For 
instance, for the case of SN1987A, this is ~ 10 46 J. 

Substituting our solution eq. (|2.10| ) to the luminosity bound in eq. ( |4.2|) we obtain 



the new astrophysical constraint on the neutrino transition magnetic moment, 

< io~V (4 3) 

~ (B/lG)J(L /lkm) x (R core /10k m y 



Thus one sees that, even for a reasonable r.m.s. field B > 10 8 G, this gives a 



more stringent bound than obtained from SN1987A [16]. It is also a more stringent 
bound than obtained above from nucleosynthesis in the presence of a primordial r.m.f. 
generated at the electroweak phase transition. 

Note that the validity of eq. Q4.3| ) requires averaging over small random r.m.f. 
distributions, so that the domain size L must be much less than the core radius 
R CO re and this, in turn, much less than the damping length Tj 1 . Indeed one can 
rewrite eq. (OF) in such a way that both of the conditions 



L < R core < ri 1 = [(8/3) ^B) 2 L ]- L (4.4) 

are verified, provided L <C R CO re- The bound shown in Fig. 2 corresponds to 
L ~ 10 meters in eq. (]4.3|). 

5. Discussion and conclusions 



We have derived bounds on transition neutrino magnetic moments connecting active 
to sterile Majorana neutrinos v a and v s for neutrinos which propagate in the pres- 
ence of random magnetic fields. We treated both the constraints that follow from 
nucleosynthesis as well as those that follow from supernova energy loss arguments. 
In both cases our constraints are typically much stronger than in the case without 
random magnetic fields. The allowed region of parameters corresponding to each case 
are shown in Fig. 1 and Fig. 2. In the nucleosynthesis case our bounds are consistent 
with the hypothesis of seeding of galactic magnetic fields by primordial fields gen- 
erated at the electroweak transition, and surviving past the re-combination epoch. 
In the supernova case, even for a reasonable r.m.s. field B > 10 8 G, our constraint 
is more stringent than obtained from SN1987A. It is also a more stringent than the 
nucleosynthesis bound obtained in this paper in the presence of a primordial r.m.f. 
generated at the electroweak phase transition. 

Our results were derived from a simple Schrodinger-type evolution equation ignor- 
ing neutrino interactions. Our results were derived for the general case of Majorana 
neutrino transition magnetic moments, but in the approximation of negligible mixing. 



They are identical to the ones derived from nucleosynthesis in ref. [jTj| for the case of 
diagonal Dirac neutrino magnetic moments. This suggests that kinetic theory effects 
neglected here are not too important and should encourage one to perform a complete 
study in which the mixing between the two Majorana neutrino species is taken into 
account. We plan to come back to this question elsewhere. 
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Figure Captions 



Fig.l. 

Constraints on active to sterile neutrino transition magnetic moments derived 
from early universe nucleosynthesis in the presence of a strong random magnetic field. 
This limit is given as a function of the domain size and field geometry parameter p, 
explained in the text. The region corresponding to the seeding hypothesis of galactic 
magnetic fields is indicated by the hatching. 



Fig.2. 



Constraints on active to sterile neutrino transition magnetic moments derived 
from supernova cooling arguments. This limit is given as a function of the random 
magnetic field strength for typical values of neutrino energy, core density and size, 
and lepton abundances. Note however that, for the dynamo mechanism suggested in 
ref. |13[] the seed field B should be larger than B > 2 x 10 8 G. 
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